Bistability of cell-matrix adhesions resulting from non-linear receptor-ligand dynamics 
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Bistability is a major mechanism for cellular decision making and usually results from positive 
feedback in biochemical control systems. Here we show theoretically that bistability between un- 
bound and bound states of adhesion clusters results from positive feedback mediated by structural 
rather than biochemical processes, namely by receptor-ligand dissociation and association dynamics 
which depend non-linearly on mechanical force and receptor-ligand separation. For small cell-matrix 
adhesions, we find rapid switching between unbound and bound states, which in the initial stages 
of adhesion allows the cell to explore its environment through many transient adhesions. 



In cell-matrix adhesion, cells use integrin-based con- 
tacts to collect information about their environment 0. 
The integrin receptors form two-dimensional clusters in 
the plasma membrane and are stabilized by connections 
to the cytoskeleton. Their extracellular domains rc- 
versibly bind ligands from the extracellular matrix like 
fibronectin, which present the ligand epitope through a 
polymeric tether. In order to efficiently explore a newly 
encountered surface and to commit itself to adhesion only 
if an appropriate combination of local signals is present, 
it is favorable for cells to start by establishing many small 
and transient adhesions sites. In the following we intro- 
duce a simple theoretical model which predicts that small 
cell-matrix adhesions are characterized by bistability be- 
tween bound and unbound states due to well-established 
principles of receptor-ligand dynamics. For biochemical 
control structures, positive feedback is well-known to re- 
sult in bistability and switch-like behavior, for example 
in the cell cycle and in the MAPK-cascade Q. Here we 
show that bistability in biological systems can also result 
from structural processes. 

Fig. H shows the model studied in the following. The 




FIG. 1: Schematic representation of our model for an adhesion 
cluster with N t bonds, out of which i are closed at a given 
time (here Nt = 5 and i = 3) . The force transducer at the top 
and the ligand tethers at the bottom are modeled as harmonic 
springs, with rest lengths It and lb and spring constants kt 
and kt, respectively. The positions of the unloaded springs 
(dashed lines) are separated by the unloaded receptor-ligand 
distance I. The extensions of the loaded transducer and bond 
springs are denoted x t and xt = I — xt, respectively. 



mechanical properties of the cell envelope holding the re- 
ceptors and of the polymeric ligands are represented by 
harmonic springs. We consider a situation in which N t 
receptor-ligand pairs are arranged in parallel along the 
cell-substrate interface. At a given time t, there are i 
closed and Nt — i open bonds. The dynamic variable i 
ranges from (completely unbound state) to Nt (com- 
pletely bound state). The probabilities for the N t + 1 
states are denoted by pi(t) (0 < i < N t ) and their time 
evolution is described by a one-step master equation: 



dp t 
dt 



r i+ ip i+ i + 9i-iPi-i - (n + 9i)pi ■ 



(1) 



This equation states that i can either decrease due to 
rupture of a closed bond (reverse rate r{) or increase due 
to binding of an open bond (forward rate gi). 

The reverse rate is = ik ff(i) because any of the i 
closed bonds can be the next to break. The single bond 
dissociation rate increases exponentially with force ac- 
cording to the Bell equation, k Q ff(i) — kQe Fb ^' l ^ Fo , with 
the unstressed dissociation rate k and the internal force 
scale F The force Fh{i) acting on each of the closed 
bonds changes dynamically and follows from mechanical 
equilibrium for the arrangements of springs depicted in 

Fig.m 
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l + i(k b /k t ) 



(2) 



Here, kh represents the stiffness of the polymeric tethers 
carrying the ligands and kt represents the effective stiff- 
ness of the cell envelope. The distance I is the average 
distance between receptors and ligands if the cluster is 
completely unbound. Eq. J5J shows that the more bonds 
are closed, the less force acts on the single bonds because 
it is shared between all closed bonds in the cluster. 

The forward rate is gi — (N t — i)k on (i) because each 
of the (N t — i) open bonds can be the next to close. The 
single bond association rate increases as the receptor- 
ligand distance Xb{i) decreases and can be written as 



k n (*) 



_ k on -i kb x b (i) 2 /k B T 

Z 



(3) 



where Z is the partition function for a ligand in a har- 
monic potential confined between — lb and Xb(i), k on is 
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the single bond on-rate if receptor and ligand are in close 
proximity, and Xb(i) = Fb(i)/kb follows from Eq. J3J|. In 
our model, temperature T essentially controls ligand mo- 
bility: the larger T, the broader the ligand distribution. 
In the limit T — > oo (high ligand mobility, k on (i) = k on ) 
and kb 3> k t (soft transducer, = k t £/i), our model 

simplifies to a case which has been used before to study 
the stability of adhesion clusters under force 0. How- 
ever, only the full model introduced here results in bista- 
bility. 

In order to demonstrate that our model leads to bista- 
bility, it is sufficient to consider the mean number of 
closed bonds, N(t) = Neglecting deviations 

from the mean, the time evolution of N(t) is described 
by the ordinary differential equation 

^L=g(N)-r(N) (4) 

where reverse and forward rate are now interpreted as 
functions of a continuous argument. The same mean field 
approach has been used before to study the stability of 
adhesion clusters under force in the limits of high ligand 
mobility and soft transducer, both for constant and 
linearly rising force 0. In Fig.0 we show the results of 
a bifurcation analysis of Eq. . In Fig. [5^, stable (solid 
lines) and unstable (dashed lines) fixed points are shown 
as function of the dimensionless unloaded receptor-ligand 
distance A = I jib and for different values of the dimen- 
sionless inverse ligand mobility (3 — kbil/ksT. For small 
A, a single, stable fixed point exists at finite AT. Here, the 
adhesion cluster is bound because force on a closed bond 
is small and the density of free ligands close to the recep- 
tors is large, therefore rupture events are rare and can 
be balanced by rebinding. For large A, there is a single, 
stable fixed point at N w 0. Here, the adhesion clus- 
ter is unbound because forces are large and ligand den- 
sity at the receptors is small, therefore rupture occurs 
frequently and cannot be balanced by rebinding. The 
transition from bound to unbound proceeds via a series 
of two saddle-node bifurcations, resulting in a window 
of bistability. This bistability results from two positive 
feedback mechanisms working in parallel. First, there is 
positive feedback for rupture: as one bond breaks, force 
on the remaining bonds is increased, thus increasing their 
dissociation rate. Second, there is positive feedback for 
binding: as one ligand binds, receptor-ligand distance 
is decreased and the binding rate for the other ligands 
increases. Fig. |2t also shows that with increasing lig- 
and mobility (decreasing (3), the bistable region shifts 
to larger A because the spatial distribution of ligands be- 
comes broader, thus stabilizing adhesion clusters through 
rebinding. In Fig. [2) we construct a stability diagram by 
plotting the A corresponding to the upper and lower bi- 
furcation points as a function of (3. Again we observe that 
A decreases with increasing f3. For the lower bifurcation 
we find the scaling A ~ f3~ x / 2 . For large (3 (small lig- 
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FIG. 2: (a) One-parameter bifurcation diagram showing sta- 
ble (solid lines) and unstable (dashed lines) fixed points of the 
mean field equation Eq. Q as function of unloaded receptor- 
ligand distance A = £/£), and for three values of inverse ligand 
mobility, /3 = k^/ksT = 0.5, 1 and 2. The other parameters 
are Nt = 10, k on /k = 1, k b /k t = 1 and k b i b /F = 0.1. (b) 
Stability diagram constructed from the location A of upper 
and lower bifurcation points shown in (a) as function of j3. 



and mobility), the width of the two-state-region is fairly 
constant, while it decreases for small (3 (large ligand mo- 
bility) as the two curves converge. 

To study the dynamics of adhesion clusters, the full 
stochastic model of Eq. has to be used. Bistabil- 
ity demonstrated in the bifurcation analysis of the mean 
field theory now corresponds to bimodality of the station- 
ary probability distribution. Again, coexistence of bound 
and unbound adhesion clusters is found for intermediate 
values of the unloaded receptor-ligand distance A. In the 
stochastic description, the system is able to switch dy- 
namically between these coexisting macrostates. Using 
Eq. , average binding and unbindingtim.es can be cal- 
culated as mean first passage times In Fig. [3] we 
plot average binding and unbinding times as function of 
the unloaded receptor-ligand distance A for cluster sizes 
Nt = 10 (solid lines) and 25 (dashed lines). A detailed 
analysis shows that their ratio effectively equals the ra- 
tio of occupancy probabilities of the stationary solution 
of the master equation. Therefore when binding and un- 
binding times are equal, unbound and bound state are 
occupied with equal probability, thus defining a stochas- 
tic transition point. The shaded areas in Fig. [3] show 
the windows of bistability derived from the mean field 



3 



/ / 


1 

f 


\ 

\ 

/ / 
' / 

/ 


/ / 
/ 

\ / 

r\ 

I * 

\ 

N 

N 

S 




binding 




1 


unbinding 










3 1 2 






3 4 

X 


5 



FIG. 3: Average binding and unbinding times in units of fco 
as function of unloaded receptor-ligand distance A = £/£{, for 
N t = 10 (solid lines) and 25 (dashed lines) at f3 = 1. The 
other parameters are as in Fig. |U The shaded areas marked 
by bars are the windows of bistability (which overlap in the 
darker region) from the mean field theory for Nt = 10 (left) 
and N t = 25 (right). 



approach. For Nt — 10, the stochastic transition point 
is below the deterministic range, while for Nt — 25, it 
falls into the deterministic range. The definition of a 
range of bistability is less appropriate in the stochastic 
case because here the existence of a bimodal distribution 
is completented by information on the transition times. 
For both deterministic and stochastic descriptions, the A- 
values for bistability increase with Nt, so that for a given 
A the bound state becomes the only stable macrostate 
when Nt grows. In the stochastic treatment, this leads 
to an super-exponcntially fast increase of unbinding time 
with Nt for fixed A, so that larger adhesions get effectively 
trapped in the bound state and only small adhesions will 
show bistable switching on an experimental time scale. 

We finally discuss how our results apply to the intcgrin- 
fibronectin bond. For the low affinity fibronectin-as/?!- 
integrin bond k ~ 0.13 Hz and F ~ 10 pN has been 
measured Q and we assume k on /k ~ 1. Fibronectin is 
a semi-flexible polymer with contour length L ~ 62 nm, 
persistence length L p ~ 0.4 nm and equilibrium length 
£b ~ 11 nm. Therefore the spring constant at small ex- 
tensions can be estimated to be kb = 3fc_eT /2L p L ~ 
Q.25pNnm _1 . From this we get k b i b /F ~ 0.27 and 
(3 ~ 7.1. In order to estimate the effective spring con- 
stant of the force transducer, we use the Hertz model for 
elastically deforming the cell with a cylinder of radius 
a. Then k t = 2aE/(l - v 2 ) ~ 0.27 pNnm -1 where we 
have used E ~ lOkPa for the Young modulus, v = 0.5 



for the Poisson ratio and a ~ 10 nm for the size of the 
adhesion cluster. Thus the ratio of force constants is 
kb/k t — 0.9 and bistability for Nt = 5 is found from 
the full stochastic model around A ~ 0.75. The actual 
unloaded cell-substrate distance is the unloaded receptor- 
ligand separation £ plus the ligand resting length lb, that 
is 20 nm. This indeed is the range of cell-substrate dis- 
tance at which mature cell-matrix adhesion starts [8j- 

In summary, our model predicts that the early phase 
of cell-matrix adhesion is characterized by rapid asso- 
ciation and dissociation of small adhesion sites which 
have attachment times in the range of several tens of 
seconds. In the presence of appropriate signals, commit- 
ment to firm adhesion can be induced through different 
mechanisms, including recruitment of additional recep- 
tors and changes in binding affinity, which increase the 
escape time from the bound macrostate. In this sense, 
the structural mechanism for bistability described here 
resembles the switch-like behavior often encountered in 
biochemical control structures. 
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